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Abstract 

Meilijson  and  Nadas  [1979]  have  obtained  stochastic  bounds  in  the 
convex  majorisation  sense  to  the  critical  path  length  of  a  project 
network  with  random  activity  durations.  In  this  paper  we  present  those 
results  in  a  more  general  framework  and,  using  similar  techniques,  obtain 
bounds  for  shortest  route,  maximal  flow  and  reliability  system  lifetime. 

Subject  classification:  #488  Bounds  for  stochastic  networks 
#672  Convex  majorisation  of  project  critical  path  length. 

#725  Stochastic  majorisation  of  reliability  system  lifetime. 

Consider  a  set  I  =  {1 .  n}  of  n  nodes,  the  base  set.  Let  1^, 

...,  1^  be  subsets  whose  union  is  I,  and  no  two  of  which  are  ordered  by 

inclusion;  {1^}  1  <  j  <  k  is  a  clutter  over  I.  The  blocking  clutter  to 

{I,}  is  a  clutter  J  such  that  I  J  *  <j>  for  all  r,  s,  and  J 

j  j.  Xf  r  s  j 

are  minimal  sets  with  this  property,  cf.  Edmonds  and  Fulkerson  [1970]. 

In  a  directed  acyclic  graph  or  in  a  two  terminal  network,  the  paths  and 
cuts  are  an  example  of  a  pair  of  blocking  clutters.  We  call  I  a  system, 


2 


and  {Ij}  1  <  j  <  k,  {J ^  }  1  <  j  <  l  the  paths  and  cuts  of  the  system.  Let 
a  weight  be  associated  with  each  node  i  of  system  I.  In  many 
combinatorial  optimization  problems  the  system  has  an  optimal  value 
function,  a  function  of  X^,...,  X^,  which  is  defined  by  the  clutter  of 
paths  or  of  cuts.  To  illustrate: 

-  Critical  path  of  a  PERT  network  (Elmaghraby  [1977]):  the  nodes 
represent  activities,  the  weights  activity  durations,  the  network  the 
precedence  constraints.  The  critical  path  length  is  the  shortest  time 
needed  to  complete  the  project,  given  by 

M  =  max  E  XA  (1) 

1<  j  <k  iel . 

J 

over  the  clutter  of  paths. 

-  Maximal  flow  (Ford  and  Fulkerson  [1962],  Lawler  [1976]):  the  nodes 
represent  pipelines,  the  weights  maximal  flow  capacities.  The  maximal 
flow  through  a  network  from  source  to  sink  is: 

L  =  min  £  X.  (2) 


l<j<Jt  ieJ, 


over  the  clutter  of  cuts. 


-  Shortest  route  (Ford  and  Fulkerson  [1962],  Lawler  [1976]):  the  nodes 
represent  sections  of  routes,  the  weights  their  lengths,  the  network 
their  connections;  the  shortest  route  from  source  to  sink  is  given  by  L, 
over  the  clutter  of  paths. 

-  Reliability  system  lifetime  (Barlow  and  Proschan  [1975]):  the  nodes 
represent  components,  the  weights  their  lifetimes.  The  system  lifetime 
can  be  expressed  in  terms  of  the  paths  {1^ }  or  the  cuts  {J  }  as: 

T  =  max  min  X^  =  min  max  X^  (3) 

l<j  <k  ielj  l<j<£  ieJ^ 

The  formulation  of  M,  L,  T  via  clutters  applies  equally  well  to 
structures  more  general  than  networks,  e.g.  precedence  relations  among 
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project  activities  can  be  defined  by  any  partial  order,  and  a  reliability 
system  can  be  defined  by  any  Boolean  coherent  structure  function.  The 
equality  (3)  holds  for  any  pair  of  blocking  clutters,  cf.  Edmonds  and 
Fulkerson  [1970]. 

The  stochastic  behaviour  of  the  optimal  value  functions  M,  L  and  T 
is  introduced  as  follows.  Let  the  weights  X^,...,  X^  be  random 


variables,  with  marginal  distribution  functions  F^,...,  F^  and  a  joint 

distribution  P.  Then  M,  L,  T  are  random  variables.  It  is  extremely 

difficult  to  obtain  the  distributions  of  M,  L,  T  -  this  is  so  even  in  the 

case  where  X^,...,  X^  are  independent,  since  different  X  ^  *  s  will  in 

general  have  nodes  in  common  and  not  be  independent.  Nor  is  it  any 

+ 

easier  to  determine  single  values  such  as  E(V),  P(V  >  y) ,  E(V  -  y)  where 
V  is  any  of  M,  L  or  T  (here  z+  =  max  (z,  o) ,  z  =  (-z)+).  Let  p  denote 
the  family  of  all  the  joint  distributions  of  X^,...,  X^  with  the  given 
marginal  distributions  F^ .  F^.  The  subject  of  this  paper  is  the 


investigation  of: 


T(x)  =  sup  E(M-x) 
n(x)  =  s^ip  E(L-x) 
a(x)  =  sup  P(T>x) 
B(x)  =  sup  P(T<x) 


We  show  how  each  of  the  functions  ¥,  n,  a,  B  can  be  calculated  as  the 
solution  to  an  appropriate  mathematical  programming  problem  which  is  in 
general  substantially  easier  than  the  calculation  of  E(M-x)+,  E(L  -  x) 
or  P(T  >  x)  for  a  particular  P  ep.  The  suprema  in  (4)  are  attained  for 
every  x,  that  is,  for  every  x  there  exists  a  joint  distribution  P  for 
which  V(x)  =  E(M  -  x)  ,  and  similar  distributions  attain  the  supremura  for 
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n,  a,  g.  The  joint  distributions  which  attain  these  suprema  can  be  chosen  to 
have  a  special  form.  Define  the  inverse  of  a  distribution  function  as 

F  L(u)  =  inf  {x  |  F  (x)  >  u} 

and  let  U  be  a  uniform  random  variable  on  (0,  1).  Then  X,  , . . . ,  X  =  F, 1 

i  n  1 

( 4>^(U) ) ,  •  •  • ,  F  ^(<|>n(U)),  where  U  is  common  to  all  the  nodes,  and  the  functions 
, . . . ,  $  are  piecewise  linear,  with  a  finite,  not  exceeding  max  (k,  Si)  , 
number  of  discontinuities.  $  ,  . . . ,  <J>n  are  obtained  explicitly  from  the 
solution  of  the  mathematical  programming  problems,  together  with  additional 
structural  information  about  the  system. 

The  functions  'V,  n,  a,  B  can  be  used  to  define  random  variables  M,  L,  T 


and  T  as  follows: 


E(M  -  x) 
E(L  -  x)' 
P(T  >  x) 
P(T  <  x) 


By  their  definition  (4,  5),  M  is  convexly  greater,  L  is  concavely  smaller, 

T  (T)  is  stochastically  greater  (smaller)  than  M,  L,  T  respectively,  for  any 
distribution  Pep. 

The  definitions  of  X  stochastically  greater  than  Y  (X  >  Y)  and  of  X 

i)  1 

convexly  (concavely)  greater  than  Y  (X  >c  Y  ,  (X  Y))  are  (cf  Stoyan  [1983]): 

X  Y  <=>  Vx  P(X  >  x)  >  P(Y  >  x)  <=>  Eh(x)  >  Eh(Y)  Vh  nondecreasing 
X  Y  <=>  Vx  E(X  -  x)+  >  E(Y  -  x)+<=>  Eh(X)  >  Eh(Y)  Vh  convex  nondecreasing 

X  >k  Y  <=■>  Vx  F.(X  -  x)~  «;  E(Y  -  x)"<=>  Eh(X)  >  Eh(Y)  Vh  concave  nondecreasing 

<=>  Vx  E(x  -  X)+  <  E(x  -  Y)+<^»  Eh(X)  <  Eh(Y)  Vh  convex  nonincreasing 

<=>  -  X  <  -  Y 

c 
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We  say  that  the  random  variables  M,  L,  T,  (T)  are  convex  upper, 

concave  lower  and  stochastic  upper  (lower)  bounds  for  M,  L,  T.  Clearly 

by  (4*  5)  they  are  sharp  bounds,  in  the  sense  that  if  for  example  Z  >  M 

c 

for  every  Pep,  then  Z  >c  M.  By  the  properties  of  >^,  >^,  Eh(M) , 

Eh(L),  Eh(T),  Eh(T)  provide  bounds  for  Eh(M) ,  Eh(L) ,  Eh(T)  for  every 
Pep,  whenever  h  has  the  appropriate  monotonicity  and  convexity 
properties;  these  bounds  are  not  necessarily  sharp,  unless  M,  L,  T,  T  are 
obtained  within  p  . 

In  general,  M,  L,  T  and  T  are  not  obtained  within  p  .  If  however  the 
system  is  series  parallel,  then  there  exist  joint  distributions  in  p  for 
which  M  =  M,  or  L  =  L  or  T  =  T,  T  =  T;  bounds  for  series  parallel  systems 
are  discussed  in  section  1,  together  with  a  discussion  of  modular 
decomposition.  In  sections  2,  3,  4  we  discuss  each  of  the  optimal  value 
functions,  M,  L  and  T  separately.  We  conclude  in  section  5  with  some 
general  remarks  on  the  type  of  bounds  presented  in  this  paper,  and  with  a 
comparison  with  other  types  of  bounds  which  appear  in  the  literature. 

The  present  work  is  based  on  a  paper  of  Meilijson  and  Nadas  [1979], 
who  derived  the  properties  of  TCx).  Some  of  the  results  on  a(x),  B(x) 
have  been  previously  obtained  by  Zemel  [1982].  A  brief  summary  of  the 
present  paper  appeared  in  Weiss  [1984].  Some  related  results  and 
extensions  appeared  in  klein  Haneveld  [1982],  and  Meilijson  [1984]. 


1.  Bounds  for  Series  Parallel  Systems. 

The  pure  series  system  with  nodes  1,...,  n  has  a  single  path  1^  = 

I  =  {!,...,  n}  and  n  singleton  cuts,  =  {1},...,  =  {n}.  The  pure 

parallel  system  has  paths  1^  =  1^=  {n},  and  a  single  cut  = 

n 

(1,...,  n).  For  the  pure  series  system,  M  =  I  X  ,  L  =  T  =  min  X  . 

i=l  l<i<n  i 

n 

the  pure  parallel  system  M  =  T  =  max  X  ,  L  =  I  X  . 

1 <i *n  1  i=l  1 


For 


The  following  three  special  joint  distributions  of  X^  , . . . ,  are 
essential  in  this  paper;  they  provide  the  bounds  for  the  pure  series  and 
the  pure  parallel  systems  (U  is  a  uniform  random  variable  on  (0,  1), 

F  (x)  =  1  -  F(x)  =  P(X  >  x) ): 

The  "perfect  tracking”  distribution  P*:  X^ .  X^  =  F^  (U),..., 

Fn'1(U)' 

-  The  "max  antithetic"  distribution  P**  defined  inductively  for  n  =  2 
by  X^,  X^  =  F^  (U),  F,,  1  (1-U)  and,  given  P**  for  X1 ,  . . . ,  X^  and  Y  ^ 

=  max  X.,  X  and  Y  .  are  distributed  by  Y  ,,  X  =  F„  1  (U),  F  ^  (1-U). 

l<i<n-l  1  n  n_1  n_1  n  Vl  n 

-  The  "min  antithetic"  distribution  P***  defined  similarly  to  P**  with 

Z  ,  =  min  X,  replacing  Y  , . 

"-1  KKri  1  ”-1 

It  is  easy  to  check  that  P*  (P**)  stochastically  minimises  (maximises) 

max  X  ,  and  P*(P***)  stochastically  maximises  (minimises)  min  X  over  p 

1  <i <n  1 

by  achieving  equality  in: 

n 

max  F  (x)  <  P  (min  X  <;  x)  <  min  (1,  Z  F  (x)) 

1 <i <n  1 <i <n  1  i=l  1 

n 

max  F  (x)  <  P  (max  X  >  x)  <  min  (1,  Z  F  (x)). 
l<;i  <n  1  <i  <n  1  i=l  1 

Also,  P*  convexly  maximizes  and  at  the  same  time  concavely  minimizes 


i’l 

P  e  p: 


seen 

by 

the  following 

argument : 

For 

every  x 

n 

+  n 

n 

(  l 
i=l 

Xi 

-  x)  <  (  Z 
i=l 

vt  -  x)  + 

Z 

i=l 

<xi  -  V 

n 

_i» 

"  . 

n 

(x  - 

Z 

X.)  <  (x  - 

E  v.)  + 

Z 

(v,  -  X,) 
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On  the  other  hand,  note  that  Z  F  (u)  is  left  continuous  non- 

i=l 

n  -1 

decreasing  in  u,  so  for  given  x  we  can  choose  uq  such  that  Z  (uq)  ^  x 

n  -1 

<  Z  Fi  (uq+),  and  we  can  then  choose  for  i  =  1 . n  for  which 

f”1  (u  )  c  v,  <  FT1  (u  +) ,  so  that  Z  v  =  x.  Using  those  v  ,  for  X.  = 
ioiio  isl1  1  1 

F  (u) ,  i  =  1,...,  n  as  in  P* ,  the  above  inequalities  hold  as  equalities. 

These  properties  of  P* ,  P**,  P***  ensure  that  the  various  bounds  are 

obtained  within  p  in  the  pure  series  and  in  the  pure  parallel  case. 

_  —  -k 

Theorem  1.1:  For  the  pure  series  system,  M  and  T  are  obtained  by  P  ,  L 
and  T  are  obtained  by  P***,  and  L  <s^,  L  for  all  Pep.  For  the  pure 
parallel  system,  M  and  T  are  obtained  by  P**,  L  and  T  by  P* ,  and  M  M 

for  all  P  e  p. 

A  useful  concept  in  the  theory  of  networks  or  clutters  is  that  of 
decomposition  into  modules  (or  autonomous  sets),  as  discussed  by  Barlow 

and  Proschen  [1975]  and  by  Mohring  and  Radermacher  [1984].  Consider  a 

*  *  *  *  * 

set  1,1  cl  and  let  I,,...,  I  be  all  the  different  subsets  of  I  of 
—  1  m 

*  * 

the  form  I  n  1^,  1  <  j  <  k.  Then,  I  is  a  module  of  the  system  I  if: 

★  * 

(i)  I-,...,  I  form  a  clutter* 

I  m 

★  * 

(fi)  Whenever  1^  n  I  =  I .  it  follows  that  for  every  r,  1  <,  r  *  ra 

*  * 

there  exists  an  s,  1  <  s  *  k.  such  that  (I.  -  I  .)  u  I  =  I  . 

i  j  r  s 

★ 

The  module  I  is  called  nontrivial  if  it  has  more  than  1  and  less  than 

* 

n  nodes.  The  quotient  system  I/I  is  formed  by  replacing  all  the  nodes 

* 

of  I  in  I  by  a  single  new  node  o,  with  a  similar  replacement  in  each 
path  of  the  clutter  {I.}.  It  is  maybe  more  intuitive  to  think  of  a 
system,  module,  and  quotient  system  in  the  reverse  order:  Start  with  the 
quotient  system  and  the  module  (those  can  be  any  two  systems),  choose  a 
node  in  the  quotient  system  (node  o  can  be  any  node)  and  replace  this 


node  by  the  base  set  of  the  module;  then  augment  the  clutter  of  the 

quotient  system,  by  replacing  each  path  which  contains  o  with  m  new  paths 

*  *  * 

in  which  o  is  replaced  by  I, .....  I  .  For  a  nontrivial  module  I  ,  call  I 

1  m 

*  *  *  * 

a  modular  composition  of  I  ,  I/I  and  call  I  ,  I/I  a  modular 
decomposition  of  I. 

Mohring  and  Radermacher  [1984]  discuss  the  preservation  of  M,  L,  T 

under  modular  composition.  Let  V  represent  any  of  the  optimal  value 

* 

functions  M,  L  or  T.  For  weights  x^,...,  x^,  let  V,  V  be  the  optimal 

★ 

values  for  the  system  I  and  the  module  I  .  Then  V  can  also  be  calculated 

*  *  * 

in  steps:  Obtain  V  ,  assign  the  value  V  as  the  weight  of  node  o  in  I/I  , 

★ 

calculate  the  optimal  value  for  I/I  .  For  X^ .  random  with  joint 

* 

distribution  P  c  p  ,  V  and  V  are  random  variables.  The  distribution  of  V 

* 

can  be  calculated  in  steps:  Obtain  the  joint  distribution  of  V  ,  assigned 

rt 

to  node  o,  joint  with  the  weights  of  the  other  nodes  of  I/I  ,  and  obtain  the 

* 

distribution  of  the  optimal  value  of  I/I  for  that  joint  distribution. 

In  the  following  sections  we  prove  that  modular  composition  also 
preserves  the  bounds  M,  L,  T,  T.  We  show  for  each  of  the  optimal  value 
functions  that: 

* 

Theorem  1.2:  If  module  I  is  replaced  by  the  single  node  o,  with 

weight  Xq  that  has  as  its  marginal  distribution  the  distribution  of  the 
*  * 

bound  for  I  ,  then  the  bounds  for  I/I  and  for  I  are  identical. 

A  general  series  parallel  system  is  defined  (inductively  in  the 

number  of  nodes  n)  as  a  system  which  is  either  pure  series  or  pure 

*  * 

parallel  or  has  a  nontrivial  module  I  and  quotient  system  I/I  both  of 
which  are  series  parallel.  Combining  theorems  1.1  and  1.2  we  have: 

Theorem  1.3:  For  a  series  parallel  system  the  bounds  M,  L,  T,  T  are 
obtained  by  joint  distributions  within  p  . 

Proof:  Combining  theorems  1.2  and  1.1  provides  a  direct  construc¬ 
tion  of  the  joint  distributions  for  which  rt,  L  and  T  are  extremal. 


9 


2.  Convex  Upper  Bounds  for  Critical  Path  Length 

In  this  section  we  discuss  the  optimal  value  function 

M  =  max  X  X 

lcUk  iel . 

3 

where  I  =  {1,...,  n},  1^,...,  I  is  a  clutter  over  I,  and  .  X^ 

have  marginal  distributions  F^,...,  and  a  joint  distribution  Pep 
(the  dependence  of  M  on  P  is  suppressed  to  simplify  notation).  We  start 
by  quoting  the  results  of  Meilijson  and  Nadas  [1979]. 

Let  ^Cx)  be  defined  by: 

'f(x)  =  inf  {(  max  E  v  -  x)  +  +  E  E(X  -  v )  +  }  (6) 

v  l<j<k  iel  .  iel  1  1 

3 

and  let  Xq  =  inf  {x  |  T'(x)  >  -1}.  It  turns  out  that  the  calculation 
of  (6)  for  x  >  xq  is  equivalent  to  the  solution  of  the  following 
mathematical  program  with  a  separable  convex  objective  function  and 
linear  constraints: 

'f(x)  =  min  E  E(X .  -  v  )+  (7) 

v  i  el  1  1 

s.t.  E  v  <  x  j  =  1,  ...,  k. 

i  el  .  1 

3 

Denote  by  , . . . ,  the  Lagrange  multipliers  (dual  variables)  of  the 
constraints . 

Theorem  2.1: 

variable  M  such  that  for  all  x, 

M  M  for  all  P  e  p. 

(iii)  For  every  x  there  exists  a  P  e  p  for  which  E(M  -  x)+  =  f(x). 

(iv)  A  particular  Pep  satisfying  (iii)  is  of  the  form  X^  , . . . ,  X^  = 

Fn^(  ^(U)) ,  where  U  ~  U(0,  1),  and  (|>^  have  at 
most  k  discontinuities  and  are  linear  inbetween. 


(i)  't'(x)  =  sup  E(M  -  x) 

(li)  There  exists  a  random 


y(x)  =  E(M  -  x)  ,  and 


10 


(v)  For  every  x  and  P  as  in  (iv),  the  Lagrange  multipliers  of  (7) 
satisfy: 

X.  <  P(M  >  x,  E  X  =  M) 

3  ielj 

With  equality  if  all  ' s  are  non  atomic  (absolutely  continuous). 

(vi)  The  constant  max  E  E  (X  )  is  convexly  smaller  than  M  for  all 

1  <  j  <k  i  el . 

J 

P  e  p  ;  in  particular  it  is  <  E(M) . 

Outline  of  the  proof:  (i)  For  every  x  and  every  vector  v,  for  every 

joint  distribution  P  e  p  and  every  realisation  X^,...,  X^  drawn  from  P: 

+  +  n  + 

(M  -  x)  <  (max  E  v  -  x)  +  E  (X.  -  v.)  (8) 

j  iel  i=l 

which  shows  that  the  right  hand  side  of  (6)  is  >  E(M  -  x)+  for  every 
Pep.  Equality  to  the  supremum  over  p  follows  from  (iii). 

(ii)  Examination  of  (6)  shows  that  'i'(x)  is  convex  nonincreasing  in  x 
with  slopes  tending  to  -1  and  0  as  x  tends  to  and  °°.  Hence  ^(x) 


defines  a  random  variable  M  according  to  (5),  and  M  M  for  all  Pep. 

For  x  >  x  ,  v  which  minimises  (6)  satisfies  max  E  v  =  x,  and  so  it 

1  <  j  <k  i  el  . 

J 

minimises  (7),  and  (6)  and  (7)  are  equivalent.  The  solution  of  (7)  at 

x  ,  say  v°,  minimises  (6)  for  all  x  <  x  . 
o  o 

(iii,  iv)  Consider  the  Lagrangean  of  (7): 

n  +  k 

¥(v,  X,  x)  =  E  E(X  -  v  )  +  E  X.  (  E  v.  -  x)  (9) 

i=l  j=l  3  iel.  1 

J 

with  Xj  >  0.  The  Kuhn  Tucker  saddle  point  conditions  for  it  are: 

E  v  <  x  and  E  v  <  x  implies  X  =  0 

iel.  iel.  1  3 

J  J 

P(Xi  >  vL)  a  E  X  <  P(X  >  v  ) 
jltslj 

_  k 

P(M  >  x)  <  E  X.  <  P(M  >  x) . 

j=l  3 


(10) 


11 


$ 


l*. 


V, 

V 


n 

.v 


For  a  given  x  >  x^ ,  let  v^,  A^  , . . . ,  ^  be  an  optimal 

k 

solution  and  a  set  of  multipliers  of  (7,9).  Let  A,  .  =  1  -  £  A., 

*  1=1  J 


I.  =  <f*  a,  =  £  A  ,  i  =  1,...,  n.  The  joint  distribution  Pel) 

Krr  1  1  *  1  r 

j|iel. 


stated  in  (iv)  is  defined  by  the  functions  <J>^ ,  i  =  1,...,  n  which  for 

m- 1 


m 


m  =  1,...,  k  +  1  and  £  A.  <  u  <  £  A.  have  the  value: 

j-i  J  j-i  j 


m-1 


<t>±(  u)  = 


(1  -  ai)  +  at  (u  -  Aj)/^ 


i  e  I. 


m 


m-1 

(1-a  )  (u  -  £  A  )/ A 

l  .  ,  J  m 

1=1  J 


(ID 


i  k  I 


m-1  m 

(v)  For  A  4°*  with  probability  A  ,  £  A.  <  U  <  £  A,  in  which 

m  1  m  i  * 


m 


i=i 


1=1 


j 


-1 


case  1  -  <  <j>^  (U)  <  1,  and  by  (10)  X^  =  F^  ( (U ) )  >  v^,  for  all  i  e 


1^,  while  at  the  same  time,  <  v^  for  all  i  £  1^.  By  (8),  we  see  that 


in  this  case  M  -  x  =  £  (X  -  v  )  and  £  v  =  x,  so  M  =  £  X  >  x. 

i  el  1  1  iel  1  lei  1 

mm  m 


The  required  inequality  follows,  and  equality  for  nonatomic  distributions 
follows  similarly.  Finally,  (vi)  holds  by  Jensen's  inequality. 

Corollary  2.2.  Modular  decomposition:  Theorem  1.2  holds  for  the  function 
M. 


*  *  *  * 
Proof:  Let  I  cl  with  clutter  I,,...,  I  be  a  module  of  I,  and  let  I/I 

-  —  1  m 


be  the  quotient  system,  with  set  of  nodes  I  =  (I  -  I  )  u  {o }  and  clutter 

of  paths  1^,...,  1^.  Let  4*,  4^,  4^  and  M,  Mq,  denote  the  bounds  for 
*  * 

the  systems  I,  1  ,  I/I  respectively.  We  look  at  the  program  (7)  and  the 
two  additional  programs: 


*0( y)  =  min  £  A  E(X1  -  u^ 
u  i  el 


s.t . 


S  *  u  <  y  j  =  1, 
i  el  j 


■,  ™ 


(12) 


.  j,  *- 


'  J-V- 


^  uA 
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and : 

Vx)  =  min  E  *  E(X  -  w  )+  +  V  (w  ) 
w  iel-I 

s.t.  E_  w  <  x  j  =  !»•••,  I  (13) 

lei . 

J 

Since  ¥  (w  )  =  E(M  -  w  )+  =  E(X  -  w  )+,  'f.(x)  Is  the  bound  for  the  module 
o  o  o  o  o  o  ’  1  7 

* 

I/I  .  We  need  to  show  that  Hf(x)  =  ^(x)  for  all  x. 

(1)  'f(x)  >  ^(x):  Let  v  be  an  optimal  solution  of  (7). 

Define : 


w  =  max  E  v 

°  l<j«an  iel*  1 
J 


Wi  Vi 


1  e  I  -  I  . 


Because  I  Is  a  module,  and  v  is  feasible  for  (7),  w  is  feasible  for  (13). 

The  value  of  the  objective  function  (13)  for  w  is  E  E(X  -  v  )+  + 

*  1  1 

iel-I 

* 

fQ( wq)  .  But  {v^}  i  e  I  is  feasible  for  (12)  with  y  =  wq,  and  so  ¥  (wq) 

<  E  *  E(X  -  v  )+,  so  the  value  of  the  objective  of  (13)  for  w  is  <;  ’F(x), 
i  el 

and  therefore  ^(x)  <,  'i'(x)  . 

(ii)  f^(x)  >  l'(x):  Let  w  be  an  optimal  solution  of  (13).  Let  u  be  an 

k 

optimal  solution  of  (12),  with  y  =  wq.  Let  v^  =>  u^,  i  e  I  ,  and 

•k  k 

v^  =  w^,  i  e  I  -  I  .  Because  I  is  a  module,  v  is  feasible  for  (7).  The 
objective  value  of  (7)  for  v  is 

Z  E(X{  -  v  )+  =  E  E(X  -  w  )+  +  E  E(X  -  u  )+ 

iel11  *  1  1  *  1  1 

i el-I  i  el 


(Xt  "  +  ’fQ(wo)  =  ^(x), 


iel-I 


thus  'l'(x)  <  4^(x). 

Monotonicity : 

Corollary  2 .3:  If  X^  are  replaced  by  so  that  Z ^  >  X^,  1  =  !»•••» 


13 


_  i 

then  the  bounds  M, ,  M  obtained  for  Z, .  Z  and  X,,...,  X  satisfy  M,  > 

1*  I  n  1*  *  n  '  1  c  i 

1 

m.  ; 

Proof :  Let  ^(x),  v^j^...,  v^^  be  the  solution  of  (7)  with  Z^  replacing  I 

~  n  (13  + 

Xt.  By  Z±  >c  X  ,  <F(x)  =  Z  E(X1  -  ')  <  Y  (x) .  Minimising  (7)  with 

i=l 

X^,  we  get  'f(x)  <  ^(x)  <;  ^(x),  so  for  all  x,  E(M  -  x)  +  <  E(M^  -  x)+.  j 

Computational  Aspects:  Nadas  [1979]  discusses  the  computational  aspects  of 
solving  the  mathematical  program  (7),  which  with  its  linear  constraints  and 
separable  convex  objective  function  is  relaively  easy.  If  E(X^  -  v^)+  is 
approximated  from  above  by  £^(v^)  piecewise  linear  and  convex,  the  program 
can  be  solved  as  a  linear  program,  and  provide  an  upper  bound  for  'f(x). 

The  approximation  is  equivalent  to  replacing  each  by  an  approximating 
discrete  distribution,  and  it  can  be  chosen  so  that  0  <  £,^(v)  -  E(X^  -  v)+ 

<  6  for  any  given  6  >  o,  uniformly  for  all  v. 

In  the  project  planning  application,  the  nodes  represent  activities 
and  the  clutter  1^,...,  1^  is  defined  by  the  partial  ordering  of 
activities,  and  consists  of  all  the  paths  from  the  start  to  the  finish  of 
the  job.  In  that  case  the  program  (7)  has  the  following  deterministic 
interpretation:  Find  activity  durations  v^  , . . . ,  v^  so  as  to  complete  the 
whole  project  by  time  x  at  minimal  cost,  where  doing  activity  i  in  duration 
Vj,  costs  E(X^  -  v  )+.  This  is  the  project  cost  curve  problem,  solved  by 
Fulkerson  [1961].  The  solution  is  effected,  parametrically  for  all  x,  by 
formulating  the  dual  problem  which  is  a  minimal  cost  flow  problem,  and 
solving  it  parametrically  for  all  flow  values;  this  can  be  done  by  the  very 
efficient  out  of  kilter  method,  cf.  Lawler  [1976].  The  minimal  cost  flow 
problem  that  arises  from  the  dual  to  (7)  is:  For  any  total  flow  value  A, 
find  flows  through  the  nodes  i,  i  =  1,...,  n  which  yield  total  flow  A, 
at  minimal  cost,  that  is: 
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K' 

ii 


& 


i 


V* 

s 


7 

w 

✓ 


(i 


n 


min  Eh.  (  a. ) 
“  i=l  1  1 


(U) 


s.t.  min  E  a.  ■  A 
l<j  i  £j j 


i  =  1 1  •  •  •  |  n 


Where  is  the  blocking  clutter  of  cuts,  and  where: 


h^(  a)  3  fQ  max  (F^(t)  -  a,  0)  dt 


- 


o  F^  (u)  du.  (15) 


The  total  flow  value  A,  and  the  flows  through  the  nodes  a^,  which  are 


obtained  from  the  solution  of  (14),  are  related  to  the  X  ' s  in  (9),  (10) 
k 


through  A3  E  X.,  a,  =  E  X..  The  corresponding  values  of  x  and  the 

J-l  J  1  j | i £l j  J 


v  's  in  (7)  can  be  obtained  from  (10). 


Redesign  of  a  PERT  network:  It  is  quite  usual  when  designing  a  project 
with  a  PERT  network  to  have  a  target  date  x  for  the  completion  of  the 
project,  and  a  nondecreasing  convex  penalty  function  C(y)  for  values  M  *  y 
>  x.  For  such  a  penalty  function, 

5(C)  3  C'(x)  T(x)  +/V(y)  V(y)dy  (16) 

where  C' ,  C"  are  the  1st  and  2nd  derivatives  of  C,  is  an  upper  bound  on  the 
expected  penalty. 

For  the  target  date  x,  the  expected  tardiness  E(M  -  x)+  is  bounded 
sharply  by  'P( x) ,  and  the  solution  of  (7)  provides  a  construction  for  the 
worst  case  distribution  with  respect  to  that  tardiness.  It  also  provides  a 
host  of  additional  information  on  that  worst  case  distribution  which  can  be 
used  to  redesign  the  project.  Let  v  =  v(x)  =  v^(x) , . . . , v^( x)  be  the  values 
of  the  solution  of  (7),  X^(x),...,X^  (x)  the  Lagrange  multipliers,  and 


a±(x) 


J  lei 


Xj  (x)*  1  "  l»  *’  n* 


j 


id 


•-.vj 


‘i.YH 


■V 


«■ 


V- 


vs" 


v-: 


m 


»-'•  ^  y  aA 
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The  values  v^(x)  provide  target  durations  for  the  activities  with 

respect  to  the  general  target  date  x.  If  we  let  ¥^(x)  =  E(X^  -  v^(x))  , 

then  ¥  (x)  is  the  expected  contribution  of  node  i  (activity  i)  to  the 

* 

total  tardiness.  Similarly,  for  a  module  I  we  get  by  solving  (13)  for 
★ 

I/I  and  due  date  x,  a  value  w  (x)  which  is  the  target  duration  of  the 

o 

* 

module  I  with  respect  to  the  general  target  date  x,  and  we  can  get  ¥  A(x) 

I 

=  E(M  -  w  (x))  (obtained  by  solving  (12),  with  y  =  w  (x)),  as  the 
o  o  o 

i k 

expected  contribution  of  module  I  to  the  total  tardiness.  If  ¥^(x)  or 

¥  .  (x)  is  inserted  in  (16)  instead  of  ¥,  we  obtain  E,(C)  and  E  .  (C)  which 
I  1  I 

★ 

are  the  worst  case  bounds  on  the  expected  contribution  of  i  or  I  to  the 
penalty.  Thus  the  v^(x)  and  wq(x)  provide  a  way  of  assigning  tardiness  and 
penalties  to  each  activity  or  module  (on  the  basis  of  a  worst  case 
analysis) . 

The  values  A^(x)  provide,  for  the  worst  case  distribution,  the 
probability  that  tardiness  beyond  x  occurs,  and  that  the  longest  path  is  1^ 
(at  least  if  all  X^'s  are  continuous  random  variables),  as  stated  in 
theorem  2.1.  It  is  also  easy  to  see  from  the  proof  of  theorem  2.1  that 
a^(x)  is  the  probability  that  tardiness  beyond  x  occurs  and  that  node  i  is 
on  the  longest  path. 

* 

Similar  quantities  can  be  calculated  for  a  module  I  .  Solution  of 

* 

(12)  with  y  =  wq(x)  provides  X' s  and  a' s  within  I  .  Solution  of  (13)  for 

it 

I/I  ,  provides  by  the  value  a  ( x)  the  probability  that  tardiness  beyond  x 

* 

occurs  and  the  longest  path  passes  through  I  . 

3.  Concave  Lower  Bounds  for  Maximal  Flow  and  Shortest  Route 


where  I  =  {1,...,  n},  , . . . ,  is  a  clutter  over  I,  and  X^ , . . . ,  have 

marginal  distributions  F^,....Fn  and  a  joint  distribution  Pep.  When 
are  the  clutter  of  paths  in  a  network,  L  is  the  shortest 

route;  when  are  the  clutter  of  cuts  in  a  network,  L  is  the 

maximal  flow.  The  results  about  L  exactly  mirror  the  results  about  M  in 

section  2.  This  is  due  to  the  duality  between  the  various  pairs  of 

concepts  occuring  here:  path-cuts,  series-parallel,  min-max,  convex- 
concave,  P(X  <  x)  -  P(X  <  x) ,  and  E(X  -  x)+  -  E(X  -  x)  . 

The  function  n(x)  in  (4)  is  given  by: 

+  n  + 

n(x)  =  inf  {(x  -  min  E  v, )  +  E  E(v,  -  X . )  }  (17) 

v  l<j  <£  i  eJ  ^  1  i=l  1  i 

J 

and  for  x  <  xq  =  sup{x|  n' (x)  <  1},  (17)  is  equivalent  to 

n(x)  =  min  E  E( v  -  X  )+  (18) 

v  i  el 

s.t.  E  v  >x  j  «  1, . . . ,  £ 

ieJ. 

J 

with  Lagrange  multipliers  X^,...,  X^. 

Theorem  3.1: 

(  i)  n(x)  =  sup  E(x  -  L)+ 

(ii)  There  exists  a  random  variable  L  such  that  for  all  x,  n(x) 

=  E(x  -  L)  +  ,  and  1  £  L  f°r  all  Pep* 

(iii)  For  every  x  there  exists  a  P  e  p  for  which  E(x  -  L)+  =  n(x). 

(iv)  A  particular  Pep  satisfying  (iii)  is  of  the  form:  X^,...,  X^ 

=  ^  (U)  )»•**»  $n(U) ) ,  where  U  "•  U(0,  l),  and  have  at 

most  £  discontinuities  and  are  linear  in  between. 

(v)  For  every  x  <  xq  and  P  as  in  ( iv) ,  the  Lagrange  multipliers  of 
(18)  satisfy: 

A  <  P(L  <  x,  E  X  =  L) 

1eJj 

with  equality  if  all  F^'s  are  non  atomic  (absolutely  continuous). 


(vi)  The  constant  min  E  E(X  )  is  concavely  larger  than  L  for  all 

l<j<Jt  ieJj 

P  e  P  ;  in  particular  it  is  >  E(L)  . 

Proof:  This  is  a  corollary  of  theorem  2.1,  if  the  problem  is  reformulated 

in  terms  of  -X  ,  with  -L  =  max  E(-X  ). 

1  l<j<£  i  el j 

The  modular  decomposition  theorem  1.2  and  monotonicity  (with  respect  to 

>  )  hold  for  L,  in  analogy  with  M. 
tc 

Computational  Aspects:  The  program  (18)  has  a  separable  convex 
nondecreasing  objective  function  and  linear  constraints,  and  can  be 
approximated  by  a  linear  program,  like  (7). 

For  the  shortest  route  application,  the  solution  of  (18)  can  be 
obtained  by  using  =  -  X^ ,  and  solving  (7). 

For  the  maximal  flow  application,  when  J^,  ...,  are  cuts,  problem 
(18)  has  the  following  deterministic  interpretation:  Find  flows  (or 
capacities)  for  nodes  1=1,  n,  so  as  to  obtain  a  flow  (maximal 

flow)  of  x,  at  minimal  cost,  where  the  cost  of  flow  v^  in  node  i  is  given 

4- 

by  E(v^  -  X^ )  ,  which  is  convex  nondecreasing  in  v^  .  This  deterministic 
problem  is  very  similar  to  the  dual  problem  for  the  project  planning 
application,  given  by  (14).  It  can  be  solved  parametrically  for  all  flow 
values  x,  using  the  out  of  kilter  method,  cf  Lawler  [1976]. 

In  applications  to  shortest  route  problems  one  may  have  a  design  value 
x  and  a  convex  decreasing  reward  function  C(y)  for  values  of  L  =  y.  In 
applications  to  maximal  flow  problems  one  may  have  a  target  flow  x  and  a 
convex  decreasing  penalty  function  C(y)  for  value  of  L  =  y  <  x.  n(x)  and  L 
provide  upper  bounds  for  the  expected  shortfall  below  x,  E(x  -  L)+,  and  of 
E(C(L)).  The  solution  of  (18)  provides  similar  information  for  redesign  as 
in  the  critical  path  applications. 
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4.  Stochastic  Upper  and  Lower  Bounds  for  Reliability  System  Lifetime 

In  this  section  we  discuss  the  optimal  value  functions 

T'  3  max  min  X,  (19) 

l<j<k  iel  1 


T"  =  min  max  X,  (20) 

l<j<£  ieJj  1 

where  I  3  {l,...,n},  1^,...,  1^  and  .  are  two  clutters  over  I,  and 

X^,...,  X^  have  marginal  distributions  F^,...,  with  joint  distribution 
function  Pep.  If  1^,...,  1^  are  the  paths  and  J^,...,  are  the  cuts 
of  a  reliability  system  (defined  through  a  network  or  through  a  general 
Boolean  coherent  structure  function  as  in  Barlow  and  Proschan  [1975]),  and 
also  for  any  other  pair  of  blocking  clutters  as  shown  by  Edmonds  and 
Fulkerson  [1970],  T'  =  T" .  For  a  reliability  system,  if  nodes  1,...,  n 
represent  components,  and  X^,...,  X^  are  the  component  lifetimes  then  T  3 
T'  =  T"  is  the  system  lifetime. 

We  will  show  that  the  supremum  functions  a(x)  and  8(x)  of  (4)  are 

given  by  solution  of  the  following  linear  programming  problems: 
k 

a(x)  =  max  E  A. 

A  j=l  3 

s.t.  Z  A.  <  F.(x)  i  3  1, . . . ,n  (21) 

k 

Z  A.  <  1 

j=l  J 

Aj  >° 


8(x)  3  max  Z  y, 

y  j=i  3 

s.t.  E  (jj<F^(x)  i  3  1 .  n 

j  !  i  e  J  . 


E  \i,  <  1 

j3l  J  „ 

>  o 
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where  F^(x)  =  1  -  F^(x)  =  P(X^  >  x) .  The  analogy  with  the  programs  (7) 

and  (18)  is  seen  in  the  dual  programs  to  (21),  (22): 
n 

a(x)  =*  min  Z  F.^  (x)  v^  +  w 
v  i=l 

s.t.  Z  v  +  w  >  1  j  =  1,...,  k  (23) 

J 


w,  >  0 
n 

3(x)  =  min  Z  F.(x)u  +  w 
u  i=l  L 

s.t.  Z  u,  +  w  >  1  j  =  1,...,  I  (24) 

ieJ  . 

J 

w,u^  >  0 


The  following  theorem  is  implied  in  parts  by  Zemel  [1982] 
Theorem  4.1: 

(i)  a(x)  =  sup  P(T'  >  x),  e(x)  =  sup  P(T"  <  x) . 


(ii)  There  exist  random  variables  T  and  T  such  that  for  all  x,  a(x)  = 

P(T  >  x),  3(x)  =>  P(T  <  x),  i.e.  T  >gT  T*  and  T  cgT  T"  for  all  Pep. 

(iii)  For  every  x  there  exist  P',  P'1  e  P  for  which  P(T'  >  x)  =  o(x), 

P(T ' '  <  x)  =  8(x)  . 

(iv)  In  particular  P',  P"  e  p  satisfying  (iii)  exist  which  are  of  the 

form  X^...,  Xn  =>  F"1  (  ^(U)  )  , . . . ,  F~1(  ^(U)  )  ,  where  U  ~  U(0,  1)  and 

^  have  at  most  max(n,  k)  (raax(n,  £))  discontinuities  and  are  linear 
inbetween . 

(v)  For  every  x  and  P' ,  P1'  as  in  (iv),  the  solutions  (21),  (22)  satisfy 


P(T’  >  x,  min  x^  =  T') 

iElj 

p.  <  P  (T”  <  x,  max  x^  =  T") 

i£jj 

with  equality  if  all  F^ ' s  are  non  atomic  (absolutely  continuous). 


Proof :  (ii)  From  the  form  of  the  objective  functions  of  (21),  (22)  u(x)  is 
nonincreasing  and  0(x)  is  nondecreasing.  From  (23),  (24)  0  c  ct(x)  <  1  and 
0  <  0  (x)  <  1.  From  (23)  (24)  one  obtains  0(-®)  =  a(  °°)  =0,  and  from 
(21),  (22)  ot(-°°)  =  0(  ®)  =  1.  a(x)  and  0(x)  depend  continuously  on  F  (x)  , 
which  are  continuous  from  the  right,  hence  a(x)  and  0(x)  are  continuous 
from  the  right.  Thus  1  -  a(x)  and  0(x)  are  distribution  functions, 
defining  T  and  T. 

(iii) (iv),  We  shall  describe  the  construction  of  members  of  P  as  stated  in 

(iv) .  For  given  x,  let  X  ( p)  be  the  optimal  basic  solution  of  (24)  ((22)). 


Let  Vi  -  1 

and  let : 


k 

E 

j=l 


U£+l  1 


I  U-, 

J-l  3 


"k+1 


=  J 


Jtt-1 


a. 

l 


E  X  <  F  (x) 


0  =  E  p .  <  F  (x) , 

J 


and  define  for  each  i  =  1,...,  n  the  following  functions  for  0  <  u  cl: 
m-1  m 

If  E  X.  <  u  <  E  X.,  and  X  =f  0,  (lcmck+l): 

j-i  J  j-i  J 


m-1 


(J>*  (u)  =  { 


(1  -  a^)  +  otj  (u  -  E^  I 


E  I. 


m-1 


(1-a)  (u-  E  X)'X 

1  j-i  J  " 


Ui 


m-1  m 

and  if  E  p,  <  u  C  E  p.,  and  p  f  0,  (1  <m  U  +  1): 
j=l  3  j-l  3  ” 


♦i'(«)  =, 


m-1 

8i  (U  "  j=Zi  WJ)/wm 


i  e  J 


m-1 


0.  +  (1  -  0.)  (u  -  E  p.)/p„ 


i  k  J 
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It  is  clear  from  the  construction  that  F  (4>|(U))  and  F  ^  (^'(U)) 


are 

both  distr 

t.buted  as 

if 

U  ~ 

U(0, 

i) 

and 

we  have 

P'  ,  P"  e  p  . 

With 

m-1 

m 

m- 

-1 

m 

probability  X^ 

<“»>-  ,£,  * 

.  < 
J 

U  < 

E 

j=l 

X. 

J 

( 

y 

E  p.  <  1 
=  1  J 

!J  <i  E  u  . ) ,  and 

j=l  J 

in  that 

case 

,  (U)  > 

1  -  v  Fi< 

x) 

(4-;f 

(U) 

< 

Fi  (x)) 

so  X 1  >  x  (Xt 

<  x)  for 

all 

i  e  I  (i 
m 

e  J  ) .  But 
m 

then  min 

Xi 

> 

X 

(max  X. 

<  x)  and  so  T* 

>  X 

i  el 

i  el 

m 

m 

(T'  ’ 

<  x) .  Thus  a(x)  <  P 

(T« 

>  x) 

and 

0< 

:x) 

<  P(T'  ' 

<  x)  for  P'  and 

P.  < 

respectively,  with  equality  actually  occuring  because  of  (i). 

( i)  We  prove  ( i)  for  a(x),  the  proof  for  6(x)  is  analogous.  Define  for 

c 

each  j  =  1,...,  k,  A  =  C>  {X.  >  x},  A.  the  complement  event  of  A.,  and 

J  i  el  .  1  J  J 

J 

define  for  each  <J>  4  K  £  (1 ,  •  •  • »  k}  B^,  =  A ^  A  *7 .  Let  =  P(BK) 

for  an  arbitrary  P  e  U.  Then  P(T  >  x)  =  P(\-2  A.)  =  P(LJ  B  )  =  E  X  .  Let 

r  J  J  K  X  ^  K 

i  be  some  node  in  I.  Then  if  i  el.  and  j  e  K,  B  e  {X  >  x},  hence 

J  K.  —  l 

<  F1  (x)  i  -  l,...,  n  (25) 

where  the  summation  is  over  all  K's  such  that  i  e  1^,  j  e  K  for  some  j. 

If  we  set  up  the  program 

max  E  (26) 

K  X 

s.t.  (25)  and  E  X^  <  1,  X^  <  0  (26) 

K 

then  the  solution  is  >  P(T'  >  x)  for  all  Pep.  The  linear  program  (26) 
has  a  variable  X^  for  each  subset  K  c  {1,...,  k}.  In  particular  this 
includes  X^,...,  X^  which  correspond  to  the  singletons. 

Consider  any  K  which  is  not  a  singleton,  say  K  £  {i ,  j}.  Then  the 
variable  X^  appears  in  every  constraint  in  which  X^  or  X^  appear.  It  is 
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then  clear  (by  examining  the  dual  problem)  that  A^  need  not  be  basic  in 
an  optimal  solution,  and  so  it  is  a  redundant  variable.  Thus  (21)  has 
the  same  solution  as  (26)  and  so  a(x)  >  P(T'  >  x) .  This  completes  the 
proof . 

Modular  decomposition  and  Monotonicity:  hold  for  a(x),  B(x)  as  for  y, 

n. 

Computational  aspects:  The  main  difficulty  in  solving  the  LP's  (21),  (22), 
is  the  number  of  variables,  one  for  each  path,  which  may  grow  exponentially 
with  the  number  of  nodes.  Consider  solution  by  the  simplex  method.  Then 
one  can  avoid  the  necessity  of  handling  such  a  large  number  of  variables, 
by  using  column  generation,  as  follows:  Let  A  be  a  basic  feasible  solution 
to  (21)  and  let  v  be  its  simplex  multipliers;  assume  £  A^  <  1  and  so  w  =  0. 
This  solution  is  optimal  if  v  is  dual  feasible,  that  is 


Iv  >l  j  =  1, . . . ,  k 

I. 

J 

to  check  that,  one  can  solve  the  shortest  route  problem  with  the  weights 


v^ .  If  the  solution  is  <  1, 
enter  the  basis. 


it  gives  a  new  path  1^  for  which  A^  should 


The  solution  for  all  values  of  x  can  again  be  done  using  parametric 

LP. 


Zemel  [1982]  has  shown  that  If  the  shortest  route  problem,  min 

l<j  <k 

Z  v  ,  can  be  solved  in  polynomial  time  then  so  can  the  LP  (21),  by 

li 

employing  an  ellipsoid  type  algorithm.  In  particular,  for  a  two  terminal 
network  the  shortest  route  problem  is  indeed  polynomial.  This  is  in 
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sharp  contrast  to  the  fact  that  the  calculation  of  the  reliability  of  a 
two  terminal  network  with  independent  components  is  NP  hard,  as  shown  by 
Rosenthal  [1975],  Ball  [1980]  and  Valiant  [1979].  This  great  difference 
in  difficulty  of  calculation  between  the  independent  case  and  the  bounds 
of  the  present  paper  is  the  most  remarkable  feature  of  these  bounds. 
Redesign  of  a  reliability  system:  The  solution  of  (21,  22,  23,  24) 
provides  similar  information  for  redesigning  the  reliability  system  as  we 
had  for  M  and  L.  Typically  one  may  have  a  target  system  lifetime  x  and  be 
interested  in  P(T  >  x) .  Else  one  may  have  a  monotone  decreasing  penalty 
function  ( y)  for  T  =  y  <  x,  or  a  monotone  increasing  reward  function 

C2( y)  for  T  =  y  >  x,  with  (C' C'^  are  the  derivatives  of  C^,  C2>: 

ECq)  =  C1(x)0(x)  +_oo/X  -  C{(y)3(y)dy 

(27) 

E(C  )  =  C  (x)ct(x)  +  r  C^(y)a(y)dy 


providing  upper  bounds  for  penalty  or  reward. 

If  0(x),  a(x)  are  replaced  by  u^(x)  F^(x)  and  by  v^(x)  F^(x)  in  (27) 
where  u(x),  v(x)  are  the  solutions  to  (24),  (23),  one  gets  an  assignment 
of  penalties  or  rewards  to  the  various  nodes  for  the  extreme  cases. 
Similar  assignments  are  obtained  for  modules.  The  probabilistic 


interpretations  of  , 

section  2. 


8^  can  also  be  used  in  redesign,  as  in 


5.  Discussion 

In  this  paper  we  have  developed  various  bounds  of  various  types  for 
the  optimal  value  functions  M,  L,  T.  It  is  fortunate  that  each  of  those 
bounds  is  exactly  of  the  type  most  useful  in  application.  The  bounds  for 
T  are  stochastic,  so  that  given  a  target  date  x  we  have  upper  and  lower 
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bounds  on  the  system  reliability  P(T  >  x)  at  the  date  x.  The  bound  for  M 
is  an  upper  bound  in  the  convex  majorisation  sense,  which  in  the  critical 
path  length  application  gives  an  upper  bound  to  the  expected  tardiness 
beyond  a  target  day  x,  E(M  -  x)+.  The  bound  on  L  is  a  lower  bound  in  the 
concave  majorisation  sense,  which  in  the  maximal  flow  application  gives  an 
upper  bound  on  the  expected  shortfall  of  the  flow  below  a  target  level  x, 
E(x  -  L)  +  . 

For  aesthetic  reasons,  or  for  some  unforeseen  applications  one  would 
nevertheless  desire  the  bounds  not  obtained  here,  i.e.  stochastic  bounds 
on  M,  L,  convex  lower  and  concave  upper  bounds  on  M  and  L  respectively, 
and  convex  and  concave  bounds  for  T.  In  the  following  discussion  we 
explain  why  these  are  unlikely  to  possess  the  same  nice  properties  as  the 
bounds  already  obtained. 

In  section  1  we  introduced  three  special  join!  distributions,  the 

* 

perfect  tracking  distribution,  P  ,  and  the  max  and  min  antithetic 

distributions,  P  ,  P  ,  which  achieve  the  various  bounds  for  the  pure 

series  and  the  pure  parallel  system.  We  state  the  following  simple  lemma: 

Lemma  5.1:  Let  Mbe  a  family  of  distributions.  If  M*  e  JH  satisfies 

M*>  M  for  all  M  e  ))1,  then  there  is  no  member  M  e  nt  for  which  M  >  M  for 
c  b  r 

all  M  e  Ht  unless  M  =  M* . 

—  —  *  —  —  *  — 

Proof:  Assume  existence  of  M.  Then  M  M  >  M,  so  M  >  M  >  so 
-  ST  c  ’  c  c 

E  (X  -  x)+  =  E^(X  -  x)+  for  all  x,  and  M  =  M* . 

Since  P*  does  not  in  general  stochastically  majorize  or  stochastically 

rainorize  I  X^  (the  only  exclusions  one  can  think  of  are  when  n-1  of  the 

variables  are  deterministic  and  p  has  only  one  member),  stochastic  bounds 
on  M  (L)  in  the  purely  series  (purely  parallel)  case  will  not  be  obtained 
within  p  . 


.  •  ••  .V  .v  A  _*• 
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We  do  not  know  how  to  construct  convex  lower  bounds  or  concave  upper 
n 

bounds  for  £  X..  Such  bounds  if  they  exist  within  will  minimize  the 
i=l  1 
n 

variance  of  £  X  and  will  provide  the  optimum  distributions  for  variance 
i=l 

reduction  in  Monte  Carlo  simulation,  as  discussed  by  Hammerseley  and 
Mauldon  [1950],  Handscomb  [1958]  and  Whitt  [1976]. 

Next  we  note  that  M  is  nondecreasing  convex,  L  is  nondecreasing 
concave  and  T',  T"  are  nondecreasing  but  neither  convex  nor  concave.  The 
generalization  of  the  bounds  from  the  pure  series  or  pure  parallel  case  to 
general  series  parallel  systems  and  the  property  of  modular  decomposition 
follow  from  the  fact  that  >  ,  >  ,  >  are  preserved  respectively  by  non- 

C  K.  ol 

decreasing  convex,  nondecreasing  concave  and  nondecreasing  functions.  We 
cannot  however  expect  the  same  to  hold  for  concave  bounds  on  M,  convex 
bounds  on  L  or  convex  and  concave  bounds  on  T,  T' .  The  stochastic  bounds 
on  T  if  obtained  within  p  are  of  course  also  sharp  convex  and  concave 
bounds.  If  they  are  obtained  outside  p  we  do  not  know  how  to  construct 
sharp  convex  or  concave  bounds. 

The  bounds  discussed  in  the  present  paper  provide  one  way  of 
circumventing  the  impossible  problem  of  calculating  the  exact  distributions 
of  M,  L  or  T.  We  conclude  by  mentioning  other  approaches  that  appear  in 
the  literature.  These  are  based  on  special  models  for  which  the 
distributions  of  M,  L  and  T  are  tractable,  namely: 

-  Markovian  systems:  If  the  weights  have  independent  exponential 
distributions  then  M,  L,  T  have  phase  type  distributions  in  the  sense  of 
Neuts  [1981],  and  their  distributions  can  be  evaluated,  though  at 
considerable  computational  effort,  e.g.  Kulkarni  and  Adlakha  [1984]. 


-  Series  Parallel  systems  with  independent  weights:  The 
distributions  of  M,  L,  T  are  obtained  by  convolutions  products  and 
complementations . 

"Perfect  tracking"  (Nadas  [1979]):  Each  weight  X^  is  of  the  form 

=  a^  +  b^  Z,  i  *  1,...,  n,  where  Z  is  a  single  random  variable,  common 

to  all  the  X^'s.  If  p  =  P(Z  <  z)  ,  then  the  p  percentiles  of  M,  L,  T  are 

obtained  from  the  values  x^,  =  a^  +  b^  z. 

By  approximating  F,  ,  ...»  F  with  phase  type  distributions  one  can 
1  n 

presumably  get  bounds  on  M,  L  or  T  using  the  first  approach. 

Series  parallel  systems  form  the  basis  for  calculations  of  bounds  on 
values  of  P(V  >  y)  where  V  is  M,  L  or  T  -  see  Robillard  &  Trahan  [1977], 
Shogan  [1977],  Devroye  [1979]  for  PERT  and  shortest  route  applications,  and 
Barlow  and  Proschan  [1975],  Esary  &  Proschan  [1970],  Shogan  [1978],  Natvlg 
[1980]  for  reliability  applications. 

The  approach  of  the  present  paper  utilises  the  third  approach  of 
perfect  tracking. 
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